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ABSTRACT 


This  document  is  the  Final  Report  on  ARF  Project  1158A01,  sponsored 
by  the  Bureau  of  Yards  and  Docks.  The  report  covers  the  second  year  of 
the  program. 

The  objective  of  the  program  is  an  analytical  and  experimental  investi¬ 
gation  of  nuclear  radiation  streaming  in  ducts  and  shelter  entranceways. 

In  the  program  just  completed  the  energy  dependence  of  duct  attenuation 
was  extended  for  gamma  rays.  Further  duct  geometry  dependence  was  also 
measured.  Thermal  neutron  attenuation  was  determined  for  a  non-point  source 
and  the  neutron  attenuation  was  found  to  be  extremely  small. 

Use  of  albedo  theory  to  successfully  describe  the  radiation  streaming  of 
gamma  rays  was  further  strengthened  and  a  computer  code  was  prepared  to 
describe  the  streaming  in  a  straight  duct.  The  code  will  represent  a  tremen¬ 
dous  labor  saving  device  for  future  calculations.  An  analytical  analysis  of  a 
straight  duct,  using  albedo  theory,  is  included  in  the  appendix  along  with  an 
outline  of  future  research  needs. 
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RADIATION  streaming  IN  SHELTER  ENTRANCEWAYS 


I.  INTRODUCTION 

In  recent  years  increasing  emphasis  has  been  placed  on  survival 
oriented  research  projects.  Both  the  Military  and  the  Civil  Defense  are 
concerned  with  the  protection  of  personnel  and  property  against  nuclear 
attack  In  terms  of  hardened  bases  and  personnel  shelters.  This  project 
deals  with  the  nuclear  radiation  aspect  of  the  weapon's  output  and  how  this 
radiation  behaves  when  allowed  to  fall  on  a  shelter  entraneeway  or  duct. 

^Vhea  nuclear  radiation  is  incident  on  a  semi-inAnlte  piano  surface* 
containing  the  entrance  to  a  duct  er  shelter  entraneeway*  a  certain  quantity 
of  the  radiation  wlU  be  transmitted.  Tho  quantity  of  the  radiation  trans¬ 
mitted  will  depend  On  tho  duet  cross  sectional  area*  number  of  bends*  the 
bend  angle*  the  length  of  each  leg*  the  wall  material,  the  energy  and  type 
of  radiation  and  the  nature  of  source*  1.  e. ,  plane  parallel*  plane  isotropic 
or  point. 

Investigation  of  the  above  parameters  is  the  objective  of  this  program. 
This  document  is  a  final  report  On  the  work  done  under  Contract  NBy-318f* 
sponsored  by  the  Bureau  of  Yards  and  Pocks*  and  summarizes  the  results 
of  the  second  year  scope  of  the  contract  covering  tho  period  October*  I960 
to  Tune*  1961. 

Since  the  inception  of  the  program  in  October.  1959.  major  emphasis 
has  been  placed  on  obtaming  experimental  data  under  sufficiemtly  realistic 
armour  BEStARCH  SOUNDATIOM  0»  IlltMOIS  IWjtltOTS  Of  tfCMIlOlOOt 
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conditions  so  that  direct  shelter  and  duct  design  criteria  might  be  formu¬ 
lated.  Existence  of  such  measurements  then  allows  analytical  investigation 
to  determine  a  general  recipe  to  obtain  the  attenuation  offer  by  any  given 
design  of  duct  or  entranceway.  The  need  for  such  research  has  been 
repeatedly  stressed  by  others  and  need  not  be  amplified  here. 

The  Foundation  is  pleased  to  have  the  opportunity  to  conduct  this 
Investigation  for  the  Navy  and  the  individual  staff  members  on  the  project 
look  forv^td  to  continuing  programs  to  complete  the  study  in  full. 
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n.  SUMMARY  OF  PRIOR  WORK 


In  the  previous  report,  ARF  1158-12,  we  summarized  most  of  the 
analytical  and  experimental  work  done  by  others.  We  have  recently 
uncovered  an  old  report  containing  differential  neutron  albedo  information. 
This  work,  done  by  Dr,  S.  W,  Kash,  who  is  now  a  member  of  our  staff,  is 
summarized  briefly, 

A,  Differential  Neutron  Albedo  Data 

Interesting  information  about  neutron  roflection  frona  different 
tnatsrials,  which  can  easily  be  related  to  the  neutron  number  albedo,  was 
reported  by  S.  W.  Kash^«  Kash  used  Chandrasekhar’s  distribution  exprea* 
•ion  for  the  reflection  of  a  constant  net  flux  falling  perpendicularly  on  a 
semi-infinite,  isotropic  scattering,  absorbing  medium.  On  integrating 
overall  perpendicular  outward  directions,  the  axpressiens  for  the  total 
fraction  of  neutrons  reflected  from  the  semi- infinite  medium,  for  an  arbi¬ 
trary  incident  angle  of  a  parallel  beam  is 

G)»  1  -  (I  C) 

where 

Coa"^  a  angle  of  incidence 

C  ratio  of  scattering  cross  section  to  toted  cross  section 

of  the  semi-infinite  medium 

H{  /^gjG)  ■  Chandrasekhar’s  reflection  distribution  function. 


.Kash,  S.  W. ,  "Reflection  oi  Neutrons  by  a  Semi -Infinite.  Isotropic 
Scattering,  Absorbing  Medium.  "  NAA-SR-275,  Reactor  Physics 
Quarterly  Progress  Report,  May-,TttJy.  1953, 
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This,  as  well  as  the  single  scattering  fraction  expression  f^(  J^^,  C),  are 
plotted  in  Figure  1.  It  should  be  noted  that  these  expressions  lack  the 
angular  scattering  dependence  in  their  derivation.  Also,  it  is  interesting 
to  note  that  the  single  scattering  approach  is  a  closer  approximation  for 
highly  absorbing  materials. 
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t'iguro  U  fraction  of  Incident  Neutrons  which  arc  Reflected 
from  a  Seml-Inflnlte,  Isotropic  Scattering  Absorbing  Medium 
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III.  ANALYTICAL  PROGRAM 


A.  Introduction 

The  problem  of  nuclear  radiation  transport  through  ducts  and 
shelter  entranceways  may  be  approached  in  many  different  ways.  The  most 
practical  approach  is  that  which  invokes  the  albedo  recipe.  The  dose  albedo 
Is  defined  as  the  ratio  of  the  dose  emerging  from  a  surface  to  the  dose 
Impinging  upon  the  surface*  For  this  approach  to  be  satisfactory,  angle 
and  energy  dependence  of  the  albedo  is  needed*  Monte  Carlo  calculated 
incident  angle  and  energy  dependent  albedo  values  are  available  for  gamma 

2 

raya  only,  and  with  very  little  information  about  the  re-etnissloa  distribution  , 
Ho  such  information  is  available  for  neutron  albedos. 

The  albedo  approach  is  the  main  analytical  method  adopted  in 
most  of  the  calculations  done  at  ARF  and  has  been  programmed  for  the 
AHF  Unlvae  1105.  The  program  has  also  been  prepared  in  Fortran  and 
vdll  calculate  the  gamma  ray  dose  and  neixtron  fluie  distribution  along  the 
centerline  In  a  straight  rectangular  duct  of  any  half  width  to  length  raUOt 
and  for  any  source  and  detector  centerline  position* 

One  of  a  number  of  reasons  for  coding  this  caleulatloit  lies  in  the 
fact  that  the  computation  of  a  rectangular  duet  Involves  a  double  integration 
(see  Appendix  III)*  and  requires  many  man  houra  of  calculation  of  effort* 
especially  when  ^/L  is  not  small  and  no  simplo  approximation  for  evaltiating 
the  integrals  Is  realistic. 

Berger,  M*  and  D.  Raso,  "Backscattering  of  Qamma  Rays*  "  N.  B.S* 
Report  No.  5982  (July  24,  1958J* 
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B, 


Albedo  Approach 


The  computer  code  utilises  the  following  analytical  approach. 

In  a  rectangular  duct  of  length  ^  and  half  width  £>,  as  shown  in  Figure  2, 

let  the  source  and  detector  centerline  positions  be  at  a  distance  and 

lu  respectively  from  the  entrance  of  a  duct.  For  a  point  source  strength 

hS  ,  the  direct  contribution  to  the  detector  is: 
o 

kSp 

The  reflection  contribution  is  to  be  treated  for  cosine  or  isotropic  re¬ 
emission  distribution  or  a  combination  of  both. 

1.  Cosine  Re -emission  Distribution 

The  reflection  contribution  to  the  detector  from  a  dif« 
ferential  roof  area  dA  located  at  tVi 

where 

.  t  2  2„ 

m  (y  m  Lj)  i  &  sec  0 
Rj  =  (y  •  sec  0 

A  is  the  fraction  of  the  total  re-emission  distribution  which  is  cosine. 

Is  the  albedo  based  on  the  first  reflection 
dA*  e  dxdy  sSsec^  0  dOdy. 
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To  include  the  multiple  reflection,  should  be  replaced 


by: 


1  -  '  y’ 


where 


J  +4  6^  y  (/^-yJ^+45* 


The  total  reflection  contribution  from  the  whole  duet  auvrase  area  to  the 
detector  is; 


kS 

0  =  — I 

c  -2 

tT 


26A  o' 


(1  -  K^e'  )Rj 


6^eec^0  dOdy 

«l*‘2 


x=  0 


2,  Isotropic  Rc- emission  Distribution 

For  this  case  the  reflection  contribution  from  an  element 
area  dA'  at  (y,  &,  0)  is 


kS  .  ^ 

■ir  ’‘rr-R-*  "i*®- 

R,  1  ZttR 

i  it 


3 

•  Simon,  A*  and  C.  £.  Clifford,  "The  Attenuation  of  Neutrons  h/ Air 
Duets  In  Shields,  "  Nuc.  Sci.  and  Engr.  156»166  {19561* 
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and  the  total  reflected  contribution  to  tha  detaotor  Is: 


T 


2 

2 


where  B  is  the  fraction  of  the  total  re-exnlssion  distribution  which  is 
isotropic,  and  B  ■  1  -  A. 


where  only  the  first  three  terms  have  been  retained.  Here,  both  K  andK, 

0  • 

are  derived  for  cylindrical  geometry  with  the  Simon  and  Clifford  assumption 
that  the  multiple  reflection  integrand  dominates  for  the  radiation  scattered 
points  in  a  plane  perpendicular  to  at  /•  It  is  also  assumed  that  these 
functions  are  not  very  seiisitive  to  the  difference  between  rectangular  and 
cylindrical  geometry  and  so  they  are  used  with  a  radius  equal  to  the  half 
width  of  the  duct*  Xn  addition  it  is  assumed  that  they  are  also  not  sensitive 
to  the  position  of  the  source  and  tlie  detector. 
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J 


The  total  flux  seen  hy  the  detectoi*  is  the  direct  plus  the 


Ireflected  cdntributions 


kS 


kS 


o 


2SA 


(1>A) 


tl-K  o'jR,  (l-K.a'J; 

C  A  X  J 


y«0  e=o 


2  2 

a*5  sec  QdOdy 


C«  Albedo  Data 

1.  Gamma  Rays 

Idonte  Carlo  calculated  gamma  ray  albedo  values  {or 
Incident  angles  ranging  from  0*  to  90*  and  for  energies  from  0. 2  mev  to 
2  tnev  were  fitted  to  a  simple  expression.  For  concrete,  the  expression  is{ 


-bcosa  ^  Best  fit  to  gamma 

■  ae  4*  e  cos  0  and  neutron  incident  (llI-C-1] 

angle  dependence. 

where  0  is  the  Incident  an^e  and  a,  b  and  c  are  parameters  to  he  evaluated 

for  the  given  gamma  ray  or  neutron  source  energy. 

2 

From  the  referenced  report ,  little  information  about  the 
re-emlsslon  distribution  is  available.  yThere  the  cosine  and  isotropic  rc- 
emission  distributions  dominate,  for  perpendicular  and  grantng  tncidenee 
respectively,  the  ^atrlbution  is  found  to  be  an  approxlraato  fit  to  the 
following: 
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A(0l  «  Jxios  0 


tIII-C-2) 


Be  st  fit  gamma 
Stixemissipn  data 

E(0J  *  I  •  y*coB  0  cni-c-3) 

wher*)  A  and  B  represent  the  cosine  and  isotropic  re^emiBsion  distribution 
parts  respectively.  It  is  assumed  that  the  albedos  £ov  other  materials  will 
behave  in  a  similar  way  as  concrete. 


2.  Neutrons 

Little  incident  angle  dependence  data  for  neutrons  Is 
available.  Squatlon  X1I*C>1  may  be  used  for  any  desired  incident  angle 
dependence  or  If  in  equation  nX-G-t»  b  e  e  ■  0,  then  a*  e  a  and  the  albedo 
ie  constant  for  all  incident  angles. 

For  re»emls8ion  dlstributionSf  again  either  isotropic 
or  coeine  dependence  should  be  used  until  dependences  such  as  equations 
IXI«C-2and  lII-C-3  for  gamma  rays  have  been  determined. 


D.  Programming  Notes 

The  following  are  Included  to  aid  Users  of  the  computer 

program. 

1.  Xis  measured  from  the  centerline  of  the  duct. 

Z.  y,  azid  are  measured  from  the  duct  entrance. 

3.  All  source  and  detector  positions  are  on  the  centerline. 

4.  The  program  assumed  symmetry  about  the  centerline  for  the  four 
faces  and  for  the  two  halves  of  any  face, 
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5.  This  program  will  work  directly  for  a  square  cross  section  tunnel, 


and  with  a  simple  modification  for  other  cross  sections. 


6. 


T. 


No  approsdmation  based  on 


is  used. 


This  program  can  be  used  for  a  circular  cross  section  duct  of  radius 


§  by  integration  along  y  (let  0^  s  0,  ^  **  ■  4)  and  • 

by  multiplying  the  reflection  contribution  part  by  N'  *  (See  Input  data) 
E,  Code  Input  Data 

The  following  are  the  input  data  symboJjs  necessary  for  alt 
cases  that  this  program  can  handle,  expressed  in  IT  (tlnivae)  and  Fortran 
(IBM)  program  language; 


Fortran 

VV 

11 

n* 

number  of  Increments  of  0  to  be  Integrated  over 
but  uvhere  n  is  aiurays  taken  ae  an  even  number. 

14 

N: 

number  o£  increments  of  y  to  be  Integrated  over 
but  where  N  Is  always  taken  as  an  even  number. 

V202 

V201 

0  : 

Q 

the  Initial  value  of  9  eaproased  In  radians. 

V203 

Y202 

the  final  valve  of  9  expressed  In  radians. 

Y204 

Y203 

the  lnlti£d  value  of  y. 

Y205 

Y204 

the  final  value  of  y 

Y206 

Y205 

kS<, 

where  kS^  is  the  source  strength  expressed  in  R/hr. 

Y207 

V206 

9 

Ax 

the  total  length  of  the  duet. 

Y208 

Y207 

ft: 

the  half  width  of  the  duet. 

Y209 

Y208 

the  source  distance  from  the  entrance  end  of  the 
duct. 

armour 
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yortran 

IT 

• 

Y210 

Y209 

4' 

the  detector  distance  from  the  entrance  end  of 
the  duct. 

Y2U 

Y210 

a: 

a  paranaeter  for  the  albedo  expression. 

Y212 

Y211 

b: 

a  parameter  for  the  albedo  expression. 

Y215 

YE12 

c: 

a  parameter  for  the  albedo  expression. 

Y500 

C499 

S‘* 

a  parameter  designed  to  be  0|  1  and  2  for  an 
isotropic,  cosine  or  both,  distributions 
respectively. 

Other  symbols  used  in  the  XT  program  are  listed  below  in 

the  order  that  they  appear  la  the  projgram. 

V2i6 

A9 

Y224 

6* 

Y225 

• 

Y226 

Y227 

• 

Y228 

R,(«  « 

0.  y 

■  0) 

Y229 

0*  y 

-  01 

y230 

cos  0  (d  ■  0.  y  «  0}  . 

C408 

a'lQ  " 

0,  y 

=  0) 

Y231 

Y221 

n 

the  y  subscript  i,  1  ^VOO 

13 

the  initial  value  ior  i  or  J 

ARHOUa  «ESC«*CW  iOtlHDAltetl  of  IlilNOIt  INttITUlK  Of  tfiCMMOlOO* 
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Y217 

9. 

J 

IZ 

the  9  subscript  5»  j  ^700 

2  2 

C401 

6  sec  9 

C402 

^1 

C403 

cos 

C404 

n* 

C403 

(first  part  of  the  Kj  expression) 

V218 

for  0  (summation  by  Simpson  rule) 

Y202 

Ay 

C501 

(y  ■  0) 

Y222 

(at  y 

Y223 

Y232 

F,  Sample  Problem 

For  demonstration  as  well  as  checking  how  and  what  the  code 
does,  a  test  problem  was  prepared,  and  a  hand  calculation  performed  as  a 
check.  The  sample  problem  Is  to  calculate  the  gamma  dose  rate  at  a 
centerline  distance  of  6  feet  from  the  entrance)  of  a  6  by  fi  foot,  18  foot  long 
straight  concrete  duct.  A  0.  2  mev  source  of  strength  kS  ,  such  that 
kS^/r  a  1  was  chosen,  and  the  source  is  posltldncd  3  feet  from  the 
entrance.  The  alhedo  parameters  a,  b,  and  c  were  flitted  to  the  NSS  data 
and  found  to  he  0*  48,  2,  0,  and  0.  08,  respectively.  A  cosine  as  well  as  an 

ARHOUS!  RttCABCH  SeUSDAnoM  Of  •tllHO'i  itisTlTwIt  Qt  f6CHt40l,e9V 
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Isfitireple  Va^elzitssloa  ^stribution  was  expected,  in  whicH  case  g  -  Z. 

CSde  previous  SeetioU. )  Integration  ovec  7  points  (6  Incretnents,  and 
11  points  (10  Ineremonts^  Ay)  along  the  length  of  the  duct  are  ehosen.  The 
Input  data  needed  is  as  follows: 

n  ■  6  y,,,  *=  18  feet  X<.  •  6  feet 

N  Z 

Ns=  10  kS  /7r^=  1  R/hrxft^  a  a  0.48 

o 

0^  ■  0  rsdiaae  28  feet  b  a  0.  2 

0  ■  0.  7854  radiaue  & «  3  feet  e  >  0. 08 

n 

y  a  0  feet  Z<,  a  3  feet  8*2 

O  1 

The  code  ''tuieondltional"  output  v/Ul  appear  as  follows: 

0^  a  0. 107250^  R/hr  ■  0  fe«t  *  »8  feet  0^  •  0.019984  R/hr 

hSo 

a  ■  0. 48  h  a  2. 0  c  ®  0. 08  “Jlf“  *  ^ 

Values  of  many  of  the  Intormediate  steps  of  the  caleulatioa  may  be  printed 
out  if  "oondUiofial"  output  is  called  for. 

We  note  that  the  code  gave  a  detector  reading  of  0,^  = 

0.  1072  H/hr  of  which  the  reflected  contribution  [0^  «  0.01998  R/hx)  is 
approximately  19%  of  the  total. 

A  hand  calculation  was  performed  to  insure  the  correctness  of 
the  computer  program. 

The  XT  and  Fortran  programs  are  given  in  Appendix  IV.  The 
Fortran  program  was  translated  directly  from  the  IT  program  but  was 
not  run. 
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j  IV.  EXPERIMENTAL  PROGRAM 

I 

From  the  beginning  of  the  Radiation  Streaming  program  the  dominant 
philosophy  which  has  guided  the  experimental  planning  baa  been  to  design 
I  experiments  and  perform  measurements  which  will  be  amenable  to  theoretical 

treatment.  Because  of  the  simplicity  of  albedo  theory  we  have  persisted  in 

I 

the  contention  that  the  albedo  formalism  will  more  than  adequately  do  the  Job 
i  of  giving  the  designer  a  tool  to  determine  the  radiation  attenuation  of  a  duct. 

Because  of  the  Importance  to  the  military  (hardened  bases)  and  to  the 

i 

I 

'  Civil  defense  requirements  (personnel  shelters),  this  program  has  to  date, 

I  layad  major  emphasis  enlarge  duets  and  entraneeways. 

)  A.  Description  of  Experiments 

Figures  3,  4,  and  5  show  different  views  of  a  concrete  entrance' 

1 

vray  having  a  single  90*  bend.  The  width  and  height  of  the  duct  are  6  feet 
i  with  1  foot  thick  concrete  walls  except  where  additional  thickness  has  been 

I  added  to  ensure  against  radiation  leakage  between  source  and  detector. 

’  ThO' walls  are  built  from  Interlocking  blocks  to  avoid  cracks. 

I  More  details  maybe  seen  in  Figures  11  through  15  of  ARF'liSB>12. 

j  B.  Gamma  Doag  Rate  Distribution  MeaBurOmenta 

Ae  the  previous  program  had  investigated  Ce-60  (F  •  i*  25  mev) 

i 

and  Cs-  137  (0.  67  mev]  dose  rate  distributions  in  the  6  by  6  foot  duct,  we 

wished  to  extend  the  studies  to  higher  and  lower  photon  energies*  For  this 

purpose  Na-24  (2.  76  mev  and  1.  36  mev)  and  Au«]98  (0*411  mev)  sources  were 
1  used. 

ARHOUR  research  EOUNOATion  Of  ItllMOIt  iHSTIIUTC  OS  TCCHMOIOO* 
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1.  Sourcaa 


tChQ  gold  scntC9  ^ae  prepaxed  In  tliQ  form  of  a  cylinder, 
one  (1)  inch  In  dlametai*  aod  sbc  (6}  inches  tong.  Tha  gold  cylinder  wall 
thickness  was  0.01  inches  to  minimize  photon  self  absorption  and  scattering. 
The  gold  cylinder  was  slipped  into  a  thin  walled  (0-  017  inch)  aluminum  can 
of  slightly  larger  dimensions  and  aluminum  solder  sealed.  Such  permanent 
construotlon  allows  the  source  to  be  reactivated  at  any  time  in  the  future. 

A  five  hour  irradiation  of  the  19.  2  gram  gold  source  resulted  In  an  estimated 
8,  I  curies  of  activity  at  shutdown.  Study  of  the  Au-19S  decay  aeheme  indicates 
that  99  out  of  every  100  disintegrations  produces  the  desired  411  kev  photon. 

The  sodium  source  was  xxtade  as  follows:  NaF,  packed  to  a 

3 

measured  density  of  1.67  g/cm  by  a  Idgh  pressure  press,  was  contained  la 
an  aluminum  can  and  sealed.  The  active  length  of  the  NaF  was  6  inches  and 
of  7/8  Inch  diameter.  The  alutmlaum  can  contained  a  total  of  54. 1  grams  of 
Na. 

Na>24  has  a  decay  scheme  like  Co»60  In  that  a  pair  of 
cascade  gammas  result  from  each  dlstategratton  ■  2. 76  mev  and 
B,  ■  X.  37  mev).  In  order  to  be  eure  that  self  absorption  was  not  modifying 
the  epeetrum,  the  source  was  given  a  one  minute  activation  In  the  Armour 
Research  Reactor.  One  hour  was  allowed  to  be  certain  that  the  2.  3  minute 
aluminum  activity  had  decayed.  The  Na-Zd  source  was  then  placed  1  meter 
from  a  3  X  3  Inch  Nat  crystal  and  on  axis  and  the  spectrum  swept  cut  with  a 
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256  cli£^pnel  ana,lyzerv  The  spectrum  observed  v^as  Igund  tO  tk  {>a]rfect 
Na-24  spectrum. 

A  twelve  hour  irradiation  of  the  54.  I  grain  sodium  source 
Inthci  oeattftl  exposure  tube  of  the  Armour  Research  Reactor  produced 
approximately  7. 4  curies  of  Na  activity.  Fifteen  hours  of  decay  was  allowed 
80  that  known  impurities  la  the  aluminum  container  were  decayed  out.  When 
the  saeaeurements  started  the  source  strength  was  appvosdmately  4.  Z  curieit. 

The  source,  in  all  eases  to  be  reported,  is  located  in  the 
exact  geometrical  center  of  the  duet,  being  raised  by  a  pulley  system  from 
a  moveable  pig  on  the  floor  of  the  duct.  Detector  oxpocure  times  were 
measured  with  a  stop  watch. 

2,  Detectors 

As  for  all  previoas  work,  the  Landsveck  ionisation 
chambers  were  again  used  for  all  doss  measuremebts.  Three-sixteenth 
(3/16)  inch  thick  bakellte  sleeves,  provided  by  Landsverl^  were  used  for 
all  measurements .  The  sleeves  were  modified  so  as  to  covet  the  top  of 
the  chambers  as  well. 

No  energy  dependence  measureznenta  were  made.  We 
assumed  tbe  chamber's  response  to  be  linear  up  to  2,  76  mev.  Figvire  6 
is  a  photograph  of  the  chambers  employed,  including  the  neutron  detector 
and  cadmium  sleeve  discussed  later. 
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Photograph  of  Detectioa  Equipment. 


3.  Hcistults 


a.  Sodium  source  measurementa 

As  v/e  Dvishsd  to  extend  the  study  to  Include  the  effect  of 
long  leg  lengths,  the  6  by  6  foot  concrete  entrsncew^y  u^ae  modified  as 
shown  in  Figure  1,  Here,  part  of  the  P^eactor'a  biological  shield  uras  used 
to  extend  the  length  of  leg  A-  ^  ^  Is  defined  as  the  distance  from  the 


source  to  the  beginning  of  the  bend, 


.  yP  is  seen  to  be  16  feet. 

m 


Figure  7  shows  the  three  source  positions  and 

Soused.  These  source  positions  define  three  ° 

and  c  7  feet.  The  positions  labeled  A,  1,  la*,  la,  2,  3,  4,  etc. ,  up 

to  11  a  are  the  detector  positions  at  which  measurements  were  made. 

Two  points  about  the  measured  data  need  explanation. 
First,  referring  to  Figuto  7,  when  a  measurement  is  made  at,  say 
position  3,  that  position  is  charaeteriaed  as  a  measurement  for  which 
J^,  ■  8.  Wo  say  the  same  measurement  would  be  obtained  if  the  remaining 
B  feet  of  the  duct  behind  position  3  were  net  present.  This  wa»  established 
in  the  first  program  and  Is  valid  for  gamma  rays  only,-  not  neutrons. 

Second,  sodium  has  a  half  life  of  approximately  15  hours 
and  thus  decayed  substanUally  over  the  week  period  required  to  make  the 
measurements.  Decay  corrections  were  computed  for  each  of  the  130 

measurements  made  for  the  three  source  positions.  The  correction  was 
made  back  to  the  time  of  the  first  measurement  (position  A)  for  the  first 
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Figuro  7*  damma  Dosa  Rato  Plstrilnition  Mea.sttrementa 
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eource  position  u9e3«  Multiple  measurement#  were  made  at  each 

detector  position  and  averaged. 

The  data  oi  Table  I  are  plotted  in  Figures  8  and  9;  the 
data  of  Table  II»  In  Figures  10  and  11:  and  that  of  Table  111  in  Figures  12 
and  13. 

Observation  of  Figures  9,  11,  and  13  shews  an  almest 
1/R^  falloff  In  the  first  leg.  This  is  quite  reasonable  in  view  of  the  low 
albedo  for  sodium  energies.  Note  that  for  neutrons  (Figure  26}  the  albedo 
appears  to  play  a  more  slgidfleant  role  in  view  of  the  nearly  l/R  falloff. 

Of  signifieance  is  the  slope  of  the  second  leg  (  data. 

If  the  second  leg  data  are  normalized  and  plotted,  all  three  *  7,  10,  and 
14}  fall  on  top  of  each  other.  Hence,  the  rate  of  falloff  in  the  second  leg 
appears  to  be  independent  of  the  length  of  for  sodium  energies. 

In  Figures  12  and  13  the  data  point  at  27  feet  is  clearly 
a  bad  measurement.  A  review  of  the  measurement  data,  however,  offers 
no  Justifteation  for  disgarding  the  point. 

Finally,  no  error  measurements  are  indicated  as  tbs 
error  la  iass  than  the  height  of  an  error  bar  that  eould  be  drawn.  From  the 
experimental  point  of  view  the  total  error  in  a  measurement  is  probably 
less  than  5%  Only  the  energy  dependence  of  the  detectors  at  2. 76  msv  is 
In  doubt. 
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TABLE  I 


Gamma  Oose  Rate  Measuremantfl  tn  6  6  3^oot  Conc:rete  Cntvancway 

Na-24  Point  Source  at  S. . 


Measurement 

position 

Centerline 
distance 
from  source 
t£t) 

(ft) 

1 

(ft) 

Dose 

rate 

(mr/min) 

Dose’^ 

attenuation 

factor 

A 

37 

14 

as  A  SB 

0.  0067 

822 

1 

36 

14 

16 

0.  0077 

715 

la 

34 

14 

14 

0.0105 

523 

2 

32 

14 

12 

0.0138 

398 

3 

28 

14 

8 

0.  0323 

171 

4 

24 

14 

4 

0.  105 

52.  2 

4a 

23 

14 

3 

0.  164 

S3.  6 

4b 

22 

14 

2 

0.  276 

20.  0 

4e 

21 

14 

1 

0.727 

7.  58 

5 

20 

14 

5.00 

1.  10 

5a 

18.5 

14 

M  •  SS 

S.  46 

6 

17 

14 

•  •• 

5.  51 

£a 

15.  5 

14 

6.65 

7 

14 

14 

7.91 

8 

11 

14 

MW 

12.4 

9 

8 

14 

MW 

22.  2 

10 

S 

14 

— 

55.  1 

^  • 

11 

3 

14 

— 

146 

lla 

2 

14 

300 

Doae  attenuation  faetor  Ib  always  defines  as  ratio  of  Position  No.  6 
reading  to  the  position  of  measurenient  reading. 
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Flgttro  8.  Measured  Gamma  Dnse  Rate  Centerline  Distribution 
in  6  by  6  Foot  Concrete  Bntraneeway  <with  Single  Right  Angle  Bend. 
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figure  9*  Measured  Gamma  Dose  Rate  Centerline  Distribution 
In  6  b7  6  Foot  Concrete  Sntranceway  with  Single  Right  Angle  Send, 
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TABLE  II 


Gamma  Dose  Rate  Measurements  in  6  by  6  Foot  Concrete  Entranceway 

Na>Z4  Point  Source 


Measurement 

position 

Centerline 
distance 
from  source 
(ft) 

i. 

(ft) 

A 

(ft) 

Dose 

rate 

(mr/min) 

Dose 

attenuation 

factor 

1 

32 

10 

16 

0.  0152 

587 

U< 

30.  S 

10 

14.5 

0.  0175 

511 

2 

28 

10 

12 

0.  0278 

321 

3 

24 

10 

8 

0.  0607 

147 

4 

ZO 

10 

4 

0.  210 

42.6 

4s 

19 

10 

3 

0.  330 

27.  1 

4b 

18 

10 

2 

0.  621 

14.4 

4e 

17 

10 

1 

3.61 

2.48 

3 

16 

10 

MM  ■> 

8. 06 

1.  11 

3a 

14.5 

10 

8.83 

6 

13 

10 

M  « 

8.  93 

6a 

11.5 

10 

•  «  a» 

11.4 

M*  «•••  M  m 

7 

10 

10 

«  M 

14.4 

8 

7 

10 

m  m  m 

28.6 

9 

4  . 

10 

m^m 

85.1 

9a 

2 

10 

«s»  «« 

340 

10 

1 

10 

mmm 

1348 
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Figure  IQ.  Measured  Gamma  Dose  Rate  Centerline  Distribution 
in  6  by  6  Foot  Concrete  Entraneeway  with  Single  Right  Angle  Bend. 
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Figaro  11*  Measurod  Qammo  Pose  Rate  Centerline  Distribution 
in  6  by  6  Foot  Conerete  Sntritneeway  with  Single  Right  Angle  Bend. 
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TABLE  111 


Gamma  Dose  Rate  Measurements  in  6  b/  6  Foot  Concrete  Entranceway 

Na-  24  Point  Source  at 


Measurement 

position 

Centerline 
distance 
from  source 
(It) 

A 

(ft) 

(ft) 

Dose 

rate 

(mr/min) 

Dose 

attenuation 

7actor 

1 

29 

7 

16 

— 

588 

la 

27 

7 

14 

0.  0409 

359 

2 

2$ 

7 

12 

0.  0467 

314 

3 

21 

7 

8 

0. 113 

130 

4 

17 

7 

4 

0.423 

34.8 

4a 

16 

7 

3 

0.  661 

22.  3 

4b 

1$ 

7 

2 

1.42 

10.  3 

4c 

14 

7 

1 

11.8 

1.  3 

S 

13 

7 

see*  a» 

12.9 

1. 13 

5  a 

U.8 

7 

14.2 

m  m  m  •• 

6 

10 

7 

Si  •  V 

14.7 

6a 

8.  S 

7 

■s  Si  « 

19.8 

m  IS  »  •  M  M 

7 

7 

7 

o»  w 

28.3 

8 

4 

7 

83.6 

8a 

2 

7 

340 

•  mm 

9 

1 

7 

... 

1322 

ARNOut  >cs(«ftCH  »euwa«*«os  of  ittiaoit  iNttiiuTi  o*  ftCHHOieof 

)} 
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Flgufe  12.  Measured  Oamma  2>>8e  Kate  Centerline  distribution 
in  6  by  6  Foot  Concrete  Entraneeway  with  Single  Right  Angle  Bend, 
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Point  Source  ■  Na»  24 


Figure  13  Measured  Gamma  Dose  Rate  Ceatarline  Dietrtbutioa 
in  6  by  6  Foot  Concrete  Entranceway  with  Single  Right  Angle  Bend 
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b.  Cold  Source  Measurementa 


The  uae  of  Au>l9S  (0. 41 1  mev)  as  a  source  v^as  entended 
to  extend  the  doee  rate  distribution  and  dose  attenuation  factor  studies  to 
lower  energies.  The  albedo  for  gamma  rays  is  known  to  have  an  extreme 
energy  dependence. 


Measurententt  made  with  Au>198  are  given  in  Tables  IV, 

V,  and  VI,  and  the  data  are  plotted  in  Figures  14  through  19.  Again,  because 
of  the  decay  of  the  source  *  2. 7  days)  ail  data  is  corrected  for  decay 
back  to  the  time  of  the  first  measurement.  (Position  No.  1,  ■  14],  In 

Figure  14,  the  point  nt  17  feet  is  again  a  bad  data  point  but  must  be  retained. 

Study  of  Figures  15,  17,  and  19«  again  show  that  the 
dose  rats  distribution  function  in  the  first  leg  falls  off  like  Just  as  ths 
much  higher  energy  sodium  does.  Comparison  of  the  data  In  Tables  III 
and  VI  for  positions  7  and  8  shows  that  for  sodium  (Table  II!)  the  dose  rats 
Mis  off  3.  5%  less  than  while  for  gold  the  dose  rate  falls  off  7%  less 
than  l/R^.  The  difference  between  the  albedos  for  say  2  mev  (sodium 
average)  and  0.4  mev,  would  give  the  factor  of  two  difference  depending  on 
the  angle  of  Incidence  chosen. 

4.  Analysis  and  Cemparlson 

In  an  attempt  to  show  trends  as  tho  parameters,  such  as 
and  energy  vary.  Table  VIl  summarizes  the  dose  attenuation 

ratios.  Seme  pertinent  observations  are: 
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TABLE  IV 


Camma  Dose  Rate  Measarement  In  6  b/  6  foot  Concrete  'Entranceway 

Au«19&  Point  Source  at 


Measurement  Centerline  J  Dose  Dose 

position  distance  ^  I  2  rate  attenuation 

from  source  factor 

<ft)  (ft)  (ft)  (mt/mln) 


1 

36 

14 

16 

0. 00313 

301 

la 

34 

14 

14 

0. 00413 

228 

2 

32 

14 

12 

0.  00597 

157 

3 

23 

14 

8 

0.  0143 

65,  7 

4 

24 

14 

4 

0.  0453 

20.8 

4a 

23 

14 

3 

0.0646 

14.5 

4b 

22 

14 

2 

0.  101 

9.  3 

4e 

21 

14 

1 

0.  185 

5.  1 

S 

20 

14 

•  •• 

0. 990* 

•  ■  fli » la  fli 

6 

17 

H 

esesee 

0, 940f 

mum  m  **m 

7 

14 

14 

— 

1.62 

s 

11 

14 

— 

2.68 

m  mu  "•mm 

9 

8 

14 

e»sa  • 

4.  72 

10 

S 

14 

— 

10.7 

iU 

2 

14 

•  •• 

67. 1 

m  m  m  m  ••  m 

Average  of  4  readings. 
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Centerline  dietanee  irosn  source  in  feet 

Figure  14.  Measured  Gamma  Dose  Rate  Centerline  Distribution 
in  6  by  6  Foot  Concrete  Fntranoeway  ^Ith  Single  Right  Angle  Bend 


Figare  tS.  Measured  Gamma  Dobq  Rata  Centarllne  Distribution 
In  6  by  6  Foot  Concrete  Entranceway  with  Single  Right  Angle  Bend 


TABTiE  V 


Ganuna  Dose  Rate  Measurement  In  &  1!>/  6  Eoot  Concrete  Enttanceway 

Au-198  Point  Source  at  £1^^ 


Measurement 

position 

Centerline 
distance 
from  source 
(ft) 

A 

(ft) 

A 

(ft) 

Dose 

rate 

(mr/miii) 

Dose 

attenuation 

factor 

1 

32 

10 

0.00616 

303 

la 

30 

10 

14 

0.  00934 

200 

2 

28 

10 

12 

0.  0119 

157 

3 

24 

10 

fi 

0.  0277 

67.5 

4 

20 

10 

4 

0.  0866 

21.  6 

4a 

19 

10 

3 

0.  150 

12.5 

4b 

18 

10 

2 

0.  227 

8.  23 

4c 

17 

10 

1 

0.  596 

3.  14 

5 

16 

10 

1.62 

1.  15 

6 

13 

10 

— 

1.87 

7 

10 

10 

3.  04 

•  MR  V 

8 

7 

10 

m  m 

3.83 

9 

4 

10 

— 

IS.  4 

9a 

2 

10 

.. 

67.7 
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mv/mtn 


Point  Source  -  Au-198  HU  kevj 


10  feet 
6  a  3  feet 


5  10  }&  20  2S  30  35  40 

Centerline  dietAnce  from  source  in  feet 

f'igxire  i  ^  Measured  Oamma  Uose  Hate  Centerline  Distribution 
In  6  by  4  Test  Concrete  Cntraneevvsy  with  Single  Right  Angle  Oe(i<! 
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Figure  17.  Measured  Gfttnma  Dose  Rata  Centerline  Distribution 
In  6  by  6  Foot  Concrete  Sntraneaway  with  Single  Right  ^ngla  Bend 
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table  VI 


Gamma  Dose  Rate  Measurement  in  6  by  6  Foot  Concrete  Entraneeway 

AuoHS  Point  Source  at  S_ 

7 


Measurement 

position 

Centerline 
distance 
from  source 
(ft) 

ift) 

-^3 

(ft) 

Pose 

rate 

(mr/mln) 

Dose 

attenuation 

factor 

1 

29 

7 

16 

0.0123 

249 

!a 

27 

7 

14 

0.016S 

1R2 

S 

25 

7 

12 

0. 0229 

134 

3 

21 

7 

8 

0. 0537 

57. 1 

4 

17 

7 

4 

0. 175 

17.  S 

4a 

16 

7 

3 

0.  2S1 

12.  2 

4b 

15 

7 

2 

0.376 

8. 15 

4c 

14 

7 

1 

2.44 

1.  26 

5 

U 

7 

•  • 

2.61 

1.18 

6 

10 

7 

mm 

3. 07 

7 

7 

7 

— 

5.  73 

8 

4 

7 

— 

16.4 

8a 

2 

7 

— 

67.2 
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CentcTline  distance  from  source  in  feet 

figure  18.  Measured  Camme  pose  Rate  Centerline  Distribution 
in  6  1>/  6  Foot  Concirete  Entraneeway  with  Single  Right  Angle  Bend 
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Figure  19.  Measured  Gamma  Dose  Rate  Centerline  IHstributlon 
in  6  )i>y  6  Foot  Concveta  Entraneaway  with  Single  Right  Angle  Bend, 
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1.  The  attenuation  ratio  of  a  duct  decreases  as  become* 
shorter* 

2.  The  attenuation  ratio  of  a  duct  is  larger  for  higher  energies. 

3.  The  performance  of  the  second  leg  distribution  is  less 
sensitive  to  the  length  of  ^ ^  for  low  energies. 


TABLS  Vtl 

Dose  itttenuatiaa  tlaties  Versus  Energy  and  Length  of 


Energy 

Duct 

half 

width 

Iftl 

6 

A 

m 

A 

m 

Dose 

attenuation 

ratio 

Nfa-24 

3 

7 

16 

588 

Ao-198 

8 

7  . 

16 

249 

Na-24 

3 

10 

16 

£187 

Au-198 

3 

10 

16 

303 

Na-24 

3 

14 

16 

715 

Au-198 

3 

14 

16 

301 

In  Figures  20  and  21  a  comparison  of  14a-24»  Co>60  and 
Au*-I98  in  the  first  (  jP  and  sscond  i  ^  ^  shown.  These  data  are 

for  >P^  *  ^  Bnd  S  «  3  feet.  The  Co-60  data  was  taken  from  Table  V  on 
page  SS  of  the  previous  Anal  report*  ARF-tlSS-lS.  The  S7a-24  and  Att-198 
data  was  normalized  to  the  Co-bO  data  at  position  number  six  (6}  which  is 
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FigUTO  20*  ComparlBOn  of  First  Leg  Centerliae  Gamma  Dose  Rate 

Distributions  for  Au,  Co,  and  Na 
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f'lgurti  21,  Comparison  of  Second  Leg  Centerline  Gamma  Dose  Rate 

Distributions  for  Au,  Co,  and  Na 
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is  the  center  of  the  right  angle  bend? 

,  2  , 

A  i/R  line  is  drawn  in  Figure  20  for  comparison.  There 

is  no  obvious  reason  why  the  Co~60  data  show  a  greater  deviation  from 
2 

1/R  .  We  therefore  suspect  the  data»  which  was  taken  under  slightly 
different  conditions  than,  the  Na  and  Au  data. 

In  Figure  21,  for  the  second  leg,  the  identical  slopes 
are  puezling.  This  might  be  explainable  on  the  basis  that  the  contributions 
from  albedo  scattering  from  the  walls  of  the  second  leg  is  negligible  com¬ 
pared  to  the  direct  contribution  frona  the  four  scattering  surfaces  which 
make  up  the  right  angle  bend  and  which  see  both  source  and  detector. 

A  tabulation  ot  ail  existing  ARF  gamma  ray  data  for  full 
Seale  concrete  personnel  shelter  entraneeways  Is  given  in  Table  VIII.  The 
attenuation  factor  is  seen  to  deereass  rapidly  with  a  decrease  in  ^ I  and 
to  also  decrease  rapidly  as  the  energy  decreases.  Note  the  behavior  of 
the  attenuation  ratio  of  Na  and  An  as  Jl,  varies.  This  Is  not  understood. 


G«  Relative  Worth  of  the  Four  Spattering  Surfaces  which  Form 
the  Right  Angle  Bend 

Given  a  duet  with  a  fixed  cross  sectional  area,  a  fixed  ^  and 
and  hence  a  definite  attenuation  ratio  for  gamma  rays  of  energy  £,  we 
wish  to  Investigate  possible  artificial  ways  of  increasing  that  attenuation. 

As  lead  has  a  lower  albedo  than  concrete,  an  obvloua  device  would  be  to 
cover  one  or  more  of  the  auxfaees  at  the  right  angle  bend.  Questions  of 
importance  are,  what  thichnecc  of  lead  and  which  surfaces  ? 
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Dose  Attenuation  Factors  for  Concrete  Shelter  Entranceways 


O  • 


V 

u 

d 

u 

d 

o 

VH 

4) 

£4! 


d 

0 

•tt  0 

ts 

< 


®  z 
IS  -S 

rt 

e 


a 


A| 

N-  £ 


<  - 


& 

u 


iSI 

mi 

?7 

a 

#vj 

a 

M 

a 

1 

00 

> :: 

(DO 

ir> 

> :: 

i 

00 

in 

> :: 

00 

4) 

4> 

0) 

.stH. 

■s 

H 


mi 

%n 

o 

N 

e« 

«■* 

3* 

I" 

m 

t- 

*0 

r» 

00 

O' 

CO 

r- 

iri 

• 

QO 

•4! 

.i 

!<! 

at 

• 

0 

• 

•4 

0 

» 

z 

I-* 

#4 

« 

ro 

fO 

M 

n 

n 

« 

1*1 

n 

fO 

r* 

!»■ 

O 

a 


^ 

-I  M  f* 


H 


H 


'5 


rt 


H  i  in  ^ 


V 

1 

H 


1 

H 


> 

1 

H 


tn 

1^ 

eo 

■'t' 

M 

0 

0^ 

00 

00 

fA* 

to 

0 

0 

0 

r* 

tn 

m 

'I' 

(0 

i*i 

CO 

m 

fM 

<'l 

m 


9 

d 


m 


*o 

lA 

p4 

e 

e» 


fo 

m 

o 


e 

s 


r« 

o 

M 

9 

O 


fo  «a 


o 

o 


ei 

o 


o 

o 


o 

d 


*0  «0  >£ 


a  >  a  >  >  > 

41  e  «  4!  4)  O 

s  a  i  a  a 


>  a 

--  V 

a 


tn 


>  >  >  > 

V  4>  4) 

a  a  a 


> 

4; 

a 


> 

4) 

a 


«  I 
a 


^  %o  ^  cOh* 

1^  ^  fO  ^  ^ 

tiJt  eJ  ri  el  »4  ^ 


CD  ^  ^  ^ 

3  ^  ^  ^ 


O  900 


d 


S  S  X 
& 


I 

:s 


I 


o 

o 

I 

0 

O 


•o 

I 

o 


a 

%o 

I 


t~ 

IH 

I 

n 

o 


«o 

o* 

n« 

I 


SO  CO 
O'  o» 

►M  »-l 
I  I 

3  3 
-«;  < 


ARHOUR  atSlABCH  XeUNOAItON  or  tLLIMOfft  IN$TlfUT(  O*  TtCMNOlOGY 


50 


A  four  (4)  curia  Co*60  settree  va^s  positioned  in  the  6  by  6  foot 
entranceway  as  shown  in  Figure  22  with  measurements  te  he  made  at 
positions  D,  and  D,.  A  lead  sheet,  6  by  6  feet,  was  placed  on  the  floor 
surface  The  data  are  given  In  Table  IX. 

The  results  show  that  the  addition  of  1/8  Inch  of  lead  to 
(floor)  reduces  the  dose  rate  by  about  11%.  The  addition  of  a  second 
1/8  inch  lead  sheet,  giving  a  total  lead  thlchnesa  of  1/4  inch  on  the  floor, 
did  not  further  reduce  the  dose  at  D^.  The  results  are  reasoiiable  and 
are  explainable  by  observing  the  average  distance  through  the  lead  the 
incident  and  scattered  gammas  must  pass.  Va  may  also  conotuda  that  1/8 
inch  added  to  the  ceiling  would  result  In  another  11%  decrease,  due  to 
synrunetry.  The  attenuation  ratios  shown  have  little  value  because  the  lead 
is  affecting  the  dose  rate  at  also.  We  should  have  used  another  pesltton 
in  (closer  to  tlie  source)  In  order  to  have  a  fixed  reference  (unaffected 
by  the  addition  of  lead). 

The  above  procedure  was  repeated  for  area  A|.  The  results 
are  given  in  Table  X.  The  tremendous  weight  of  the  lead  (72  x  72  x  1/8 
inch  sheet)  precluded  banging  more  than  1/8  inch  on  A^.  Clearly,  1/8 
inch  of  lead  on  A^  results  in  a  greater  attenuation  than  on  A^  or  A^.  The 
difference  (n/  2%)  In  the  no  lead  values  of  (0. 0223  •  Table  UC,  and 
0, 0217  "  Table  X)  is  due  to  reproducibUtty  of  source  and  detector.  The 
data  in  the  two  tables  tnentloned  above  were  taken  on  two  different  days, 
the  setup  being  necessarily  disassumbled  in-between. 
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Aj  (aide  wali) 
A2  (back  wall) 
Aj  (floor) 

A^  VCAillng) 


Figure  22.  Experimental  Geometry  for  Determining  Contribution 
from  Scattering  Surfaces 

which  Form  Right  Angle  Bend  of  the  6  by  6  Foot  Concrete  Entranceway 


TAB1.E  IX 


Effect  of  Lead  Sheet 

on  Scattering  Surface 

'  -  -  -  -  . 

Position 

Reading 

BNG 

Dose  Rate 

Attenuation 

%  Change 

(mr/min) 

(mr/min) 

(mr/min) 

Ratio 

in 

No  Lead 

8.  85 

0. 00065 

8.85 

o 

397 

0.  0229 

0.  00065 

0.  0223 

1/8  Inch  Lead  on  Area  A^ 

8.  96 

0. 00065 

8.96 

o 

453 

11.  2% 

0.  0204 

0. 00065 

0.  0198 

1/4  Inch  Lead  on  Area  A^ 

8.  80 

0.00065 

8.80 

6 

440 

10.5% 

°1 

0.  0207 

0.  00065 

0.0200 
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TASL£  K 


Effect  of  Lead  Sheet  aa  Scattering  Surface 


Position 

R  eading 

BNG 

Dose  Rate 

Attenuation 

%  Change 

(mr/min) 

{mr/ min} 

(mr/m!n} 

Ratio 

in  D 

1 

D  /D 

6  1 

No  L>ead 

D, 

3.  45 

Q. 0006S 

8.  45 

0 

390 

Q.0223 

0. 00065 

0.0217 

1/s  Inch  Lead  on  Area  A^^ 

•8,  31 

0.00065 

8.  31 

0 

462 

17% 

0.0166 

0.00065 

0.018 
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To  sln\ula.te  %  perfect  radiation  trap  the  lead  was  removed  and 
the  wall  represented  by  area  (ceiling)  was  removed  completely.  The  dose 
rate  at  Dj  decreased  by  18%«  Comparison  with  the  data  of  pages  55  and  61  of 
the  previous  final  report  {ARF  1158-12)  show  an  18%  reduction  also.  The 
reproducibility  here  Is  gratifying. 

Finally,  areas  and  A^  were  separately  removed  to  determine 
their  worth.  A  summary  of  the  results  is  as  follows: 

Area  Reduction  in  dose  rate  due  to  removal  of  area 


22% 

10% 

18% 

18% 


Referring  again  to  Figure  22  note  that,  for  example,  area  A,,  is  not 
defined  au  as  to  represent  the  entire  area  of  the  floor  seen  by  the  detector 
at  D^«  The  shaded  areas,  seen  by  the  source  and  detector,  will  contribute 
to  the  reading  at  along  with  A^.  Thus,  the  reductions  In  dose  rate 
represent  the  specific  total  ooi^rlbotlon  that  only  tho  designated  area  (A^) 
makes  to  The  sum  of  the  four  area  contributions  does  not  therefore, 

.  add  to  the  total  reading  at  because  of  tho  contributions  due  to  the  shaded 
areas  and  any  corner  effect  contribution* 

Xn  attempting  to  draw  conclusions  based  on  the  measured  data, 
caution  must  be  exercised.  The  measuf«i:  .^sts  shew  that  each  of  the  areas 
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reduces  the  dose  rate  at  hy  tipprOxxnuitely  20%,  except  v^hich  gitfeS 
10%.  The  measueements,  however,  *ay  nothing  afaont  the  ahsolute  quantity 
(dose)  of  radiation  from  each  area. 

Consider  versus  Ay  The  incident  angle  for  Is  0*  .  The 
incident  angle  for  A  is  73*.  At  1  mev,  the  albedo  for  0*  incidence  Is  only 
1/5  of  the  albedo  at  73“ ,  Or  consider  A  ^  versus  A^,  The  angle  of  incldenee 
Is  the  same,  but  the  angles  of  emergence  to  reach  the  detector  directly 
are  In  two  different  planes.  This  difference  cannot  be  evaluated  a*  fio  real 
differential  albedo  emergence  angle  dependence  is  known  to  exist. 

Finally,  it  must  be  emphasized  that  the  reduction  in  tha  doso 
rate  due  to  removal  of  any  one  of  the  four  areas  has  an  and  S 

dependence.  As  Jl ^  and  ^incvBaae^  the  shaded  areas  of  Figure  22  will 
decrease  and  for  a  shorter  the  effect  of  any  one  scattering  area  will 
diminish.  Also,  this  experiment  was  done  with  a  long  in  order  to 

reduce  corner  effect  contributions  to  to  a  minimum  (a^  1%), 

D,  Thermal  Neutron  Flux  Distribution  Measurement  from 
Plane  Thermal  Source 

During  the  ftrst  year  of  the  program  measurements  of  epi- 
cadmium  and  sub»cadmiunlV  neutron  centerline  flux  distributions  were  made* 
(See  pages  71  and  79  of  ARF  1158-12).  Attenuation  ratios  were  found  to  he 
small  and  the  agreement  between  measurement  and  calculation  was  poor. 
The  poor  agreement  was  believed  to  be  due  to  lack  of  any  differential 
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Zieutron  albedo  Information  needed  tor  th.e  caleulations.  Also,  evldenne  was 
obtained  which  indicated  the  neutron  number  albedo  waa  much  large;  than 
gamma  albedos  and  a  single  measurement  gave  a  value  of  0.  6'75  fo;  the  ept> 
cadmium  neutron  number  albedo. 

The  second  year  neutron  effort  attempted  to  increase  the 
accuracy  of  measurements  and  to  Investigate  the  effect  of  plane  sources. 

The  sections  to  follow  report  the  results  of  these  measurements  and  describe 
the  conditions  of  the  experiments. 

1.  Source 

As  the  full  scale  6  by  6  foot  concrete  walled  shelter 
entranceway  is  built  so  that  the  face  of  the  graphite  thermal  column  of  the 
Armour  Research  Reactor  covers  the  entire  entrance  of  leg  this  graphite 

face  can  be  used  as  a  source  of  neutrons.  The  geometry  of  the  experiment 
is  shown  In  Figure  E3. 

Aa  v.'e  were  interested  in  sub-cadmium  neutrons  (we  will 
call  them  thermal),  we  wished  to  be  sure  that  the  neutrons  coming  from  the 
graphite  face  were  thermal.  A  series  of  measurements,  using  a  cadmium 
covered  BF^  then  a  bare  BF^.  were  made.  The  ratio  of  the  bare  to  cadmium 
covered  count  rate  will  be  directly  proportional  to  thermal  to  fast  neutron 
flux  ratio.  As  we  had  no  information  as  to  the  average  energy  Of  the  fast 
(epi-cadmium  neutrons)  no  energy  sensitivity  corrections  to  the  fast  count 
were  made. 
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J^iguss  Zit  Thermal  Neutron  Oistrlliutiun  Measurements 
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Measurements  of  the  thermal  to  fast  count  rates  vrere  made 
at  positions  12,  6,  and  1  (Figure  23}  and  the  ratios 

Bare  -  Cd  —  Thermal 
Cd  Fast 

formed.  These  measurements  Indicated  that  99«  75%  of  the  neutrons  leaving 
the  graphite  thermal  column  face  and  entering  the  duct  were  sub^cadmlum 
i  ^  44  ev).  Again,  no  counter  efficiency  correction  to  account  for  the 

1/v  drop  in  the  boron  cross  section  was  possible  a*  the  average  epl» 
cadmium  neutron  energy  was  unknown.  The  results  are  sufficiently  cQn» 
elusive  so  as  to  call  the  measured  dlstributiony' thermal  or  sub«cadmluiru 

2<  Detectors 

As  was  mentioned  In  previous  reports,  no  tissue 
equivalent  neutron  dose  counters  are  commercially  available,  and  as  such, 
all  our  measurements  are  of  number  flux.  For  this  purpose  a  new  and 
more  sensitive  BFj  counter  was  purchased.  It  bad  a  4-1/2  inch  active 
length  and  was  1  inch  In  diameter  with  an  aluminum  cathode  and  60  cm  of 
Hg  pressure.  This  detector  was  checked  for  its  gamrna  ray  sensitivity.  Its 
plateau  was  determined  for  neutrons  in  the  gamma  ray  field  expected  to  be 
present  during  the  measurement  and  for  the  length  of  cable  needed.  The 
voltage  used  was  3025  volts  with  a  1  mv  sensitivity  setting  on  the  dlserimbt- 
ator.  (The  detector  can  he  seen  in  Figure  6.  ) 
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Source  Distribution 


In  order  to  obtain  the  source  distribution  a  BF  covinter 
was  covered  by  a  1/8  inch  thick  cadmium  cylinder.  Only  the  end  area  of  the 
counter  was  left  bare*  The  coimter  was  positioned  so  that  its  axis  was 
perpendicular  to  the  graphite  thermal  column  face.  With  the  reactor  at  0.  1 
watt  power  level,  the  counter  was  moved  horizontally  and  vertically  across 
the  center  of  the  graphite  face  in  six  i.nch  steps.  The  results  of  these 
measurements  are  given  in  Table  XI  and  plotted  in  Figure  24,  A  cosine 
distribution  is  plotted  for  comparison  and  the  measured  distribution  can 
rather  reasonably  be  described  as  an  isotropic  cosine  source  vertically  and 
horizontally. 

4*  Thermal  Neutron  Distribution  Measurements  -  Results 

The  reactor  power  was  again  raised  to  0. 1  watt  and  thermal 
neutron  measurements  made  at  the  24  pasitions  indicated  in  Figure  23.  AU 
tneasuxements,  except  positions  B  and  C,  were  made  on  the  centerline  of  the 
enttanceway.  Note  the  concentration  of  pointe  at  the  bend  and  the  measure* 
mcnts  at  positions  labeled  A,  B,  and  C.  the  data  are  given  In  Table  XIl 
and  plotted  In  Figures  25  and  26. 

Xn  Figure  26,  tho  data  from  Figure  Zi  itf  repletted  and 
labeled  cosine  source.  In  addition,  previous  data  for  a  point  thermal 
•ource  (ARF  I1S8-12,  Table  Xll,  page  78),  Is  plotted.  The  point  source 
'  data  was  normalized  to  give  the  same  flux  at  position  six  <6),  45,  466  counts 
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table  xt 


Tkerinal  Nautroix  Source  Distribution 
Reactor  Power  0.  1  Watt 


Position 

Distance 

Position 

Distance 

along  the 

3&  axis 

,  from  the 
left  side 
(inches) 

cts/min 

along  the 
y  axis 

from  the 
lower  end 
(inches) 

Cts/mla 

1 

c 

5. 175 

12 

0 

5.  600 

2 

6 

15.919 

13 

6 

18,  583 

3 

121 

28,  ^6 

14 

12 

31.  755 

4 

18 

40.  745 

15 

1ft 

41. 997 

5 

7A 

50,  716 

16 

24 

53.867 

6 

30 

58, 700 

6 

30 

58.700 

7 

36 

53.  945 

17 

36 

54. 173 

a 

42 

45.151 

18 

42 

42.  386 

9 

48 

32. 096 

19 

48 

34.  035 

10 

54 

19.  090 

20 

54 

18,  Bll 

11 

60 

3.519 

21 

60 

5.192 
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TABI^S  xn 


I 


Measured  Thermal  Neutron  Flux  Centerline  Distribution  in  6  by  6  Foot 
Concrete  Entranceway  from  Isotropic  Cosine  Thermal  Neutron  Source 


Position 

Centerline 
distance 
from  source 
(ft) 

A 

(ft) 

A 

(ft) 

BF 

3 

(eounts/mln) 

Number 

attenuation 

ratio 

A* 

40 

17 

„ 

1.492 

30. 47 

B4« 

17 

1,602 

28.  38 

m  m 

17 

— 

1.588 

28.  63 

1 

39 

17 

16 

1.763 

25.79 

la 

38 

17 

15 

2.  080 

20. 96 

lb 

35 

17 

12 

3.425 

13.  27 

2 

32 

17 

9 

5.248 

8.  66 

i 

29 

17 

6 

8,  469 

5.  37 

4a 

26 

17 

3 

14,  340 

3.  17 

4b 

25 

17 

2 

17,953 

2.53 

4e 

24 

17 

1 

21, 888 

2. 08 

5 

23 

17 

MM 

35, 979 

1.26 

5a 

22 

17 

MM 

40.  063 

1.14 

Sb 

21 

17 

ss  m 

42.  847 

1.06 

6 

20 

17 

•  m 

45,466 

6a 

19 

17 

mm 

48, 152 

6b 

18 

17 

m  m 

52,  350 

7 

17 

17 

mm 

5S,  446 

8 

14 

17 

mm 

80,869 

9 

11 

17 

mm 

115.321 

to 

B 

17 

mm 

160,  275 

11 

5 

17 

—  •• 

249,  732 

12 

2 

17 

.. 

672,685 

13 

0 

It 

1,616,  682 

■  IS  •» 

0 

One  foot  outside  the  second  leg 
three  feet  to  the  right  on  No.  1. 
***three  feet  to  the  left  of  No.  1. 
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JTigura  25*  MeaBured  Thermal  Neutron  Flux  Centerline  Distribution 
iti  6  by  6  Foot  Concrete  Entraneeway  with  Single  Right  ^Angle  Bend 

Aniieut  RtitSftCM  ROUHOATtON  O*  IlLIHOIS  IM4T*TUtE  OR  tSCHMOlOO* 


Pigure  26.  Compariaon  of  Thermal  Neutron  Flux  Centerline  Distributions 
for  Point  and  Cosine  Isotropic  Thermal  Neutron  Sources, 


per  miiiutd.  Note  that  the  ratios  are  not  the  same  andhenee 

Influence  any  comparison. 

For  convenience,  a  1/R  line  is  drawn  in  order  to  observe 
tho  rate  of  flux  fail  off  in  the  first  leg  for  the  two  different  sources.  Recall 
that  the  gamma  dose  rate  falls  off  aft  l/R^, 

Comparing  the  slopes  of  the  cosine  and  point  source 
eurvea  with  the  l/R  line,  both  fall  off  like  1/R  In  the  first  leg  (  ^  j).  In 
the  second  leg  the  cosine  source  curve  falls  faster  than  the  point 

SOttrce.  The  attenuation  ratios  (No.  6/No,  1}  are  quite  different,  and 
disturbingly  low.  The  cosine  source  has  the  larger  attenuation  ratio  which 
is  fortunate  since  a  weapon  detonated  at  a  distance  from  the  entrance  to  the 
•belter  would  appear  more  as  a  plane  isotropic  source,  rather  than  a  point 
iourco  at  the  entrance. 

A  number  of  points  must  be  borne  in  mind  In  attempting 
to  draw  conelualone  from  the  measurements.  The  data  is  for  thermal 
neutrons,  which  from  a  delivered  dose  viewpoint  are  grossly  less  important 
than  fast  neutrons.  However,  the  shelter  entrance  might  be  located  at  a 
dtetance  from  the  hurst  point  at  which  a  high  thermal  flux  is  present.  If 
it  is  argued  that  the  scattering  surfaces  at  the  right  angle  bend  could  be 
covered  with  beral  or  cadmium,  the  attenuation  of  the  thermals  would  be 
largo.  Unfortunately,  the  horal  or  cadmium  becomes  a  large  source  of 
hard  capturo  gamma  rays.  Zt  is  not  clear  that  we  have  made  any  gain  from 
a  delivered  dose  viewpoint. 
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Second,  Cos  the  ^oint  source  the  namhec  distribution  la 
falls  off  like  l/R,  instead  o£  l/R^  as  the  ganuna  distribution  doea, 
because  of  the  Riaj;nitude  of  the  albedo  for  thermal  neutrohs  W  0.  V).  We 
believe  the  slower  fall  off  in  0  for  the  point  source  may  be  due  to  the 
closeness  of  the  source  to  the  bend  {  ^  9  feeth  thus  Increasing  the 

corner  effect  contribution,  (  =  17  feet  for  cosine  source). 

If  the  fast  neutron  dose  attenuation  is  anywhere  near  as 
Xewas  the  number  attenuation  (p.  69,  ARF  1158>12),  a  real  problem  exists. 
Creator  research  efforts  must  be  placed  on  neutron  dose  distribution  data 
and  theory. 
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SUMMARY  AND  COI^CI-USIONS 


Comjiietioik  of  this  program  has  resulted  in  thtt  {ivailaHUty  of  a 
8u'i»8tanUal  amount  of  measured  data  on  radiation  distributions  in  duets, 

Sttsh  data  is  highly  necessary  input  to  duct  and  shelter  design  and  has  «on- 
rtderahle  value  In  that  It  serves  as  a  guide  to  the  development  of  analytical 
toehnlques. 

Speeifleally*  gamma  ray  distribution  measurements  were  extended 
In  energy  and  geometrical  dependence.  Thermal  neutron  Hux  distribution 
measurements  for  a  non  poiiU  source  were  made  and  a  computer  code 
written  which  has  a  running  time  of  approximately  two  minutes  and  which 
giver  dose  distrlbutlene  in  straight  square  ducts. 

Considerable  effort  was  put  out  in  order  that  com{dete  details  of  all 
meaeuremente  he  clear  and  that  the  calculations  made  and  included  In  thle 
report  be  equally  clear  ae  to  methodology* 

The  third  year  of  the  program  will  concentrate  on  measurements  in 
ducts  of  much  smaller  cross  section  than  6  by  6  feet  and  Include  a  seeond 
right  angle  bend  or  third  leg.  The  computer  code  will  be  extended  to 
handle  duets  containing  right  englo  bends.  The  code  will  represent  a 
slgiUficant  labor  saving  device  to  project  personnel  and  will  be  available 
for  general  use.  Persons  haring  a  computation  to  be  solved  may  prepare 
the  program  for  their  computer  or  submit  the  input  to  ARt*  for  running  here* 
The  Foundation  is  pleased  to  have  the  opportunity  to  conduct  this 
program  for  the  Bureau  of  Yards  and  DoeVa.  ^Ve  look  forward  to  the  next 
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phass  of  the  program  and  'wouid  be  equally  pleased  to  have  the  opportunity 
to  discuss  other  research  projects. 
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APPENDIX  H 

FUTURE  RESEARCH  RECOMMENDATIONS 

The  itemized  list  of  suggested  research  problems  is  brief  but  will 
be  clear  to  those  knowledgeable  in  the  field.  The  list  Is  not  in  any  order 
oi  importance,  nor  is  it  complete.  We  are  prepared  to  undertake  these 
studies  and  would  welcome  the  opportunity  to  have  suggestions  as  to  their 
conduct. 

1.  Kleasurement  of  dose  rate  distribution  In  two  leg  ducts  for  fast 
neutrons  from  point  sources  to  develop  the  theory  in  order  to 
handle  analytically. 

2.  Measurement  of  fast  neutron  dose  rate  distributions  in  2  leg  duets 
for  plane  sources* 

3*  Measurement  of  thermal  neutron  flux  distributions  in  2  leg  ducts 
for  plane  sources. 

4.  Measurement  of  fast  neutron  dose  rate  distributions  in  3  leg  duet 
(two  right  angle  bends). 

5.  Experimental  determination  of  differential  (energy  and  angular 
dependence)  neutron  albedos* 

6.  Measurement  of  gamma  dose  rate  at  position  one  (1),  resulting  from 
removal  of  neutron*  by  absorbers  at  the  right  angle  bend  surfaces* 

7.  Measursmant  of  gamma  dosa  rate  dlstrlbutiona  in  2  leg  ducts  for 

plane  isotropic  sources.  (Use  one  or  two  sources  and  employ 
superposition*  ] 
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8.  Measurement  of  gamma  dose  rate  distributions  in  three  P)  leg 
duets. 

9.  Measurement  of  differential  gamma  albedos. 

1  0.  Completion  of  the  computer  cede  to  calculate  the  dose  rate 

distribution  and  the  attenuation  of  a  duct  of  any  geometry,  any 
energy  of  source,  and  for  any  number  of  bends. 

11.  Measurement  of  gamma  dose  rate  distributions  specifically  for 
fission  product  gammas  as  a  function  of  time.  A  17-233  fission 
plate  would  be  used  as  a  source. 

1  2.  Measurement  of  neutron  and  gamma  dose  rate  distributions  lit  one 
toy  one  foot  cross  section  ducts. 
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SEMl^RIOOROUS  ALBEDO  ANALYSIS 
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APPENDIX  nt 


SSMI-RIGQRQUS  ALBEOO  AWAX^YStS 

Consider  a  straight  rectangular  duct  o£  hall  thickness  S,  length  ^ 
and  having  a  wall  thickness  eg,uivalent  to  man/  mean  free  paths  for 
acattering  of  the  source  radiation.  Let  the  sourea  ha  a  point  and  located 
on  the  duct  centerline  and  In  the  plane  el  the  entrance.  The  detector  Is 
placed  on  the  centerline  and  In  the  exit  plane  ol  the  duct. 

The  geonaetry  of  the  problem  la  shown  in  Figure  The 

differential  wall  surface  area  dA  Is  on  one  of  the  four  walls  of  the  rectangular 
duet.  The  reference  coordinate  system  is  at  the  source.  dA  is  a  distance 
R,  from  the  source  and  U,  from  the  detector.  Only  the  uneollided  and 
once  reflected  radiation  is  to  be  consideted.  This  is  an  excellent  approxi¬ 
mation  lor  gamma  rays  but  poor  lor  neutrons.  The  uneollided  flux  reaching 
the  detector  is 

The  flux  striking  dA  from  the  source  is  at  an  incident  sngld  of  f  measured 
from  the  perpendicular  to  dA  and  is  seen  to  be 


4TfR: 


CoS  ^ 


The  reflected  or  re<-emls8ion  flux  from  dA  is  a  Iraetleno^ol  that 
which  fell  on  dA.  called  the  albedo  and  is  a  function  of  the  type  and 
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enei'gy  and  angle  o£  Incidence  of  the  radiation.  Thus,  the  quantity  of 
radiation  leaving  dA  and  reaching  the  detector  0  it 


cos  0 


e^(0.  E) 

tVk\ 


where  the  re>emleslen  dlatclbutloA  it  taken  as  isotropic.  Since  we  assume 
a  moAoenergetle  aoarce  of  radiation  the  albedo  is  taken  as  0^(0}  only. 

From  Figure  ZII-1,  the  following  may  be  asoertainadi 


dA  a  dxdy 
cos  0  a  6/11^ 
r[  •  S^K^iy^ 

R?  ■  6^  ♦  +  U  -  y)* 


The  total  reflected  flux  reaching  D  from  the  one  wall  surface  is 


Considering  that  there  are  four  symmetrical  walls  and  substituting  in  the 
celatlonshlps  from  Figure  llt>1. 
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The  double  integral  is  not  conveniently  handled  analytically  without  the 
introduction  of  undesirable  approximations  (such  as  84?  «  1}  and  in 
addition  we  wish  to  Investigate  the  behavior  of  the  integrand.  Such  investi¬ 
gation  recjuires  resorting  to  double  numerical  Integration,  a  sonjewhat  tedious 
process  if  any  accuracy  is  desired. 

!n  order  to  study  the  behavior  of  the  integrand, the  following  para¬ 
meters  were  chosen  as  representative  of  a  practical  duct  such  as  a  shelter 
entranceway. 

^  =  14  feet 

i  m  2  feet 
Ax  ■  1  foot 
Ay  ■  1  foot 

S  ■  0*4tl  tnev  {Au»  198] 

c  1.0  -  1  tnr/hr*  x  1  curie 
a  ■  ROg  ■  "  "  '  ~  j 

1. 2  X  iO  y/cm  -see 

Merger  and  Raso**  eoxnputed  the  energy  and  incident  angle  dependenea 
of  gamma  ray  albedos  In  oonetete.  Thelf  data  for  0.  2  and  0*5  tnev  la  plotted 
and  a  linear  ordinate  extrapolatiyil  dune  to  the  estimate  0.  4  mev  albedos.  The 
data  is  shown  in  Flgave  Ill-2i  plotted  as  a  function  of  the  cosine  of  the  incident 
angle,  with  energy  ae  a  parameter. 

^Radiological  Health  Handbook  fp.  139]. 

’''*NBS  Report  5982. 
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figure  lU-2*  Poac  Albedos  Xor  Conerete  as  a  Function  of  the  lueident  Angle. 
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The  integrand  is  evaluated  for  15  values  of  y  and  4  different  values 
of  X  as  shown  in  the  Tables.  The  behavior  of  the  integrand  is  seen  in 
Figures  111-3  and  III-4.  (The  10^  factor  is  for  convenience.  )  Recall  that  the 
source  is  at  y  =  0  and  the  detector  at  y  =  =  14  feet.  We  see  that  the  wall 

section  closest  to  the  duct  center  (x  =  0)  and  the  general  wall  area  nearest 
the  source  (y  =  0  to  3  feet)  will  contribute  heavily  to  the  integrated  value  of 
I(x,  y).  Also,  that  the  area  near  the  detector  (y  =  11  to  14  feet)  is  more 
important  than  the  area  near  the  center  of  the  duct. 

Numerical  integration  of  the  four  areas  yields 

I(x  )  =  150.88 
o 

I(x^)  =  140.74 
l(x^)  =  116.55 
I(x^)  89.85 

In  Figure  II1-4,  I(x)  is  plotted  against  x  and  a  second  integration  again 
performed  by  the  trapezoid  rule  resulting  in  I(x,  y)  =  378.  63.  We  have 
expressed  the  source  in  equation  (5)  in  terms  of  dose,  hence  the  solution 
to  five  (5)  required  dose  albedos  and  gives  the  answer  in  dose. 


I(x,y)  dxdy 


3  X  3.  7  X  10^^ 

10^  X  1.  2  X  10^ 
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a  0.  0642  mr/min. 
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14.6287  3130.52  18  56349 


The  uncollided  dose  is 


D  =  - :r  =  0.  225  mr/min- Curie 

"  4K  ,e  ^ 

The  total  dose  rate  seen  by  the  detector,  14  feet  from  a  one(l)  curie  point 
source  of  Au-198  (E  0.  41  mev)  inside  a  iS  by  6  foot  rectangular  duct  is 


9  U 


^  0.  225  +  064 


=  0.  29  mr/min  per  curie. 

The  albedo  scattered  contribution  (D  )  to  the  total  (D_)  is  seen  to 
be  approximately  22%. 

We  now  wish  to  have  an  experimental  check  as  to  the  validity  of  the 
analytical  approach*  We  may  obtain  this  by  referring  to  Table  IV.  At 
position  7  the  source  Is  14  feet  away*  The  measured  value  is  1.  62  mr/min 
from  a  source  of  unknown  curie  strength  and  we  assume  the  wall  behind 
the  detector  to  make  a  negligible  contribution.  We  have  determined  that 
22%  of  the  reading  {s  due  to  albedo  scattering*  Therefore,  1.  62  x  0.  78  « 

I.  27  mr/mln  direct  from  the  source.  For  a  point  source  of  1  curie  of 
Au-198  we  have 


5  mr/min  at  3  feet^  or  0.  23  mr/mln  at  14  feet. 


For  the  direct  value  of  1,  27  mr/min  the  source  strength,  at  the  time 


of  the  measurement,  must  have  been 


1.  27 
0.  23 


a  5.  5  curies. 


Radtolngiral  Health  Handbook  (p.  139). 


ARHOUR  DEtlARCfl  POUNOAIlOa  OF  ILtINOIS  INSTTTUTE  OE  TECtlNOtOGt 


91 


Using  the  5.  5  curie  gold  source,  we  measure  a  dose  rate  of 
1.  62  mr/min  at  14  feet  from  the  source.  Our  analysis  yields  0.  29  mr/min- 
Curies  x  5.  5  curies  =  1.60  mr/min,  some  l%low. 

The  computer  program  described  in  Section  III  and  in  Appendix  IV, 
was  given  the  problem  with  a  six  (6)  curie  source. 

The  input  statement  appears  as 

«ll 

10 

40 

-y201 

0. 

.7854 

0. 

14. 

20.6 

14. 

3. 

0. 

14. 

.425 

2.1 

.05 

»c499 

2. 

■ENDATA 

■NOMORE 


and  the  output  as 


10656466  00 
14000000  02 
42500000  00 


0  14000000  02 

30000000  01  0 

21000000  01  50000000.01 


24037515-01 

14000000  02 
20600000  02 
0 
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I 

I 

]  The  totsil  dose  rate  at  the  detector  =  0.  10658  R/hr  =  1.  77  mr/min. 

However,  for  a  5.  5  curie  source 

i ! 
i  1 

1  5  S 

D_  =  1.77x-^=  1.62inr/hr. 

I  °  ^ 

I 

1 

Measurement  -  1.62  mr/min. 

t 

i  Hand  calculation  -  1. 60  mr/min. 

i  Code  calculation  -  1.62  mr/min. 


1 

I 


i 
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APPENDIX  IV 

IT  AND  FORTRAN  COMPUTER  PROGRAMS 
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A  UNIVAC  1 105,  IT  PROGRAM  FOR  RADIATION  DISTRIBUTION 

IN  STRAIGHT  DUCTS 


1 


700 


70© 


33  t2,k,S, 


1. 

2. 

4. 

6. 

7. 

8. 

9. 


10. 

11. 

12. 

13. 

14. 

15. 

16. 

17. 

18. 

19. 

2©. 

21. 

22. 

23. 

25. 

0 

26. 

26.5 

27. 

33 

27.2 

25 

27.3 

30 

27.4 

9 

27.5 

28 

27.6 

31 

27.8 

26. 

13 

28.6 

32 

29. 

2 

30. 

31. 

32. 

33. 

34. 

35. 

36. 

37. 

7 

38. 

39. 

4©. 

8 

41. 

42. 

43. 

6 

44. 

5 

45. 

ARMOUR 

« 

READ 

Y2l6:a(Y202-Y201)/ll 

i3:“0 

y224!-Y207*Y2©7 

y225;“Y206*-Y206 

y226:-Y208*Y208 

y227:-Y209»Y209 

y228:»Q6e,(Y226+Y224)Q 

y229:«Q6e,(Y227+Y224)Q 

y230:«»Y207/Y228 

c0;«(Y210»Q2e,-Y23©»Y211q)+(Y212*Y230) 

c4o8:««CQ 

yQ : -( Y206fY207-( Y225+2 .*Y224 )/q6e , ( Y225+4.*Y224 )Q 
)  /(2.*Y207) 

c© ; »( 2  .*Y207*Q6e,  Y230C1 )/( ( 1-YG^  C©  )n229  ) 

Y0 : aY224/( Y225+4.*Y224 ) 

Y©s-( Y2©6/( 2.*q6e, ( Y225+4.*Y224)Q ) )*( l+Y©+( 2. 25*Y 
0»Y0)) 

c0 1 -C0+( 1-Q6e , y23©q) /(1-Y0»  c40O ) 

c© t  o( c0* y205*Y224*c4W )/( Y228*Y228*Y228* Y229*Y229 ) 

y231 s  «( Y205*.7854)/( ( Y209-Y20O )*( Y2©9-Y20O ) ) 

ty231ty0tc0tc408 

i3*«l 

y221:«y203 

23ii5,i3>l,i4+l, 

I3J-1 

y217:-y2©1 

g13ify2©1>©. 

g13ify221>0. 

y217;«y217+y216 

i3:-2 

3; 12, 13,1,11+1, 

I2:ii2.1 

ci2:-Y224/(Q4E,Y217a*Q4E,Y217«) 

c4©1:oCI2 

c 1 2 : -q6e , ( C 1 2+( Y221-Y208 )*  ( Y221-Y20e ) )Q 

c402!-CI2 

ci2;-Y207/CI2 

Q7iFc499ii© 

o8ifc499"1« 

Q6tFc49^'*2, 

c403:-0 

y230:-0 

05 

cli©3s'*l'> 

y230:“1. 

05 

c403;-CI2 

Cl 2  J  a( Y210*{ q2e , ( -C 1 2*Y2n )Q ) )+( Y212*C 1 2 ) 
c4©4i«CI2 


F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

r 

F 

F 

F 

F 

F 

F 

F 

F 

F 

r 

F 

F 

r 

r 

F 

F 

F 

F 

F 

F 

F 

F 

F 


COMPUTE  H 
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46.  g19ifc403=0.  f 

47 .  Y 1 2 ; -( ( ( Y205-Y221 )* ( Y206,Y221 ) )+2 .*Y224 )/q6£ , 

48.  (((Y206-y221)»(Y206-Y221))+4.*Y224)Q  f 

49 .  Y 1 2 ; =( Y2©6+( 2 .*Y207 )-( ( ( Y221*Y221 )+2 .*Y224 )/q6e , 

50.  ((Y221->^Y221)+4.*Y224)Q)-Yi2)/(2.*Y20T)  f 

c  1 2 :  o(  2  .*Y207/(1-Y  1 2#  C 1 2 )  )*Q6e,  c403/(  ( ( Y221-Y209)*(  Y 

'52.  221-Y209))+C401)Q  f 

53.  TYl2TCt2Tl2Tc4©4  F 

54.  19  Gllir.c403til.  F 

55.  yi2s-Y224/(((Y206-Y221)*(Y206-Y221))+4.»Y22U)  f 

56.  yi2:«.1+YI2+(2.25*YI2*Yi2)  f 

57.  yi2:-YI2*(Y206-Y221)/q6e,(((Y206-Y221)*(Y206 

58.  -Y22i))+4.*Y224)Q  F 

59.  c405;bYI2  f 

60!  yi2:-Y224/((Y221*Y221)+4.*Y224)  F 

61 .  Y 1 2  s »( 1 .+Y 1 2+( 2 . 25* Y I 2*Y 1 2 ) )*( Y221/Q6e , ( ( Y221*Y221 ) 

62.  +4.*Y224)Q)  f 

62.5  tyi2tc405  f 

63.  Yi2;a(Yi2-fc405)*.5  f 

64.  17  g27ifc4©3>0.  f 

66.  26  ci2!*0.  f 

67.  27  ci2taCI2+(l.-Q6E,c403Q)/(l.-YI2*C4©4)  F 

67.5  tyi2tci2  f 

68.  11  ci2:«i(CI2*Y205*C404*C4©1)/((C4©2*C4©2*C402)*(((Y221- 

69.  y209)*(Y221-Y2©9))+C401))  f 

70.  tyi2tci2tc4©3tc4©1  f 

70.5  I2S-I2+1  F 

71.  3  y217j-Y217+Y216  f 

72.  y21O:-0  f 

73.  4,12,2,2,11-2,  F 

74.  4  Y2l8:aY2l0+(2.*CI2)+4.*C( 12+1)  f 

75.  y( i5+50O):o(y216/3.)*(c0+(4.*c1)+y218+ci1)  f 

76.  Y220J-(Y20l!-Y2e3)/l4  f 

80.  c501i-y501  f 

83.  22  c( 1 5+50®) J-y( 15+500)  F 

84.  Tc( i5+5®0)ti5  f 

85.  23  y221iwY221+Y22©  f 

86.  y222:-0  f 

87.  16  24, 15,3,2, »4-l,  f 

88.  24  y222:-y222+(2.*C(i5+50©))+4.*c( 15+501)  f 

88.5  ty222  f 

89.  y223«-(y220/3.)*(c5©1+(4.*c502)+y222+C( i4+1))  f 

90.  y232j«Y223+Y231  f 

91.  14  ty232ty203ty2©4ty223  f 

92.  15  ty206ty207ty206ty2©9  f 

93.  ty228ty229ty230tc4©8  f 

94.  ty216ty220tv501tc5@1  f 

95.  21  ty210ty211ty212ty205  f 

96.  10  t(0.)  f 

97.  '3l  FF 


Note;  There  is  no  difference  between  upper  and  lower  case  letters  in  the  IT  Language. 
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AN  IBM  1620,  FORTRAN  PROGRAM  FOR  RADIATION  DISTRIBUTION 

IN  STRAIGHT  DUCTS. 


DIMENSION  Y(701),  C(701) 

1  READ,  I1.D1,I4.D4.  C(500) 

DO  40  J=202,  213,  1 
READ,  Y(J) 

40  CONTINUE 
Y(217)=(Y(203)-Y{202))/D/ 

13=0 

Y(225)=Y(208)*Y(208) 

Y(226)=Y(207)*Y(?.07) 

Y(227)=Y(209)*Y(209) 

Y(228)=Y(210)*Y(210) 

TEMP=Y(227)+Y(225) 

Y(229)*SQR(TEMP) 

TEMP=Y(228)+Y(225) 

Y(230)=SQR(TEMP) 

Y(231]j=Y(208)/Y(229) 

TEMP»i(-Y(231))*Y(212) 

C(1)=(Y(211)+EXP{TEMP))+(Y(213)*Y(231)) 

C(409)=.C(1) 

TEMP=Y(226)+4.  ♦Y(225) 

Y(I)=(Y(207)+Y(20e)-(Y{226)+2.>l‘Y(225))/SQR(TEMP))/(2.  *Y(208)) 

C{1)=(2.*Y(208)*SQR(Y{231)))/((1..Y(1)*C(1))*Y(230)) 

Y(l)rY(225)/TEMP 

Y(1)=(Y(207)/(2.*SQR(TEMP)))*(1.+Y(1)+(2.25*Y(1)*Y(1))) 

C{  1)=C(  l)+(  1, -SQR(Y(23 !)))/( 1.  .Y(  1  )*C(409)) 

C(1)=(C(1)*Y(206)*Y(225)*C(409))/(Y(229)*Y(229)*Y(229)*Y(230)*Y(230)) 

Y(232)=(Y(206)*.7854)/((Y(210)-Y(209))*(Y(210)-y(209))) 

IF(SENSE  SWITCH  1)  41,42 

41  PUNCH,  Y(232),Y(1),C(1).C(409) 

13=1 

42  Y(222)=Y(204) 

Ml  =13+1 

M2  =14 +2 

DO  23  15«M1,M2, 1 
25  13=1 

Y(218)=Y(202) 

IF(Y(202))28,28,  13 
28  IF(Y(222))3l,31,  13 

31  Y(218)=Y(2ie)+Y(217) 

13=2 

13  M4=I3+1 

M5=Il+2 

DO  3  M3=M4,M5, 1 

32  I2=M3-1 

C(I2)=Y(225)/(COS(Y(218))*COS(Y{218))) 

C(402)=C(I2) 

TEMP5C(I2)+(Y(222)-Y(209))4(Y{222)-Y(209)1 

C(I2J=SQR(TEMP) 
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43 

44 
7 


8 


6 

5 


45 


46 
19 

47 


48 

49 

26 

27 

50 

n 


51 

52 
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C(403)=C(I2) 

C(12)=y(208)/C(I2) 

IF(C(500))43,7,43 

IF(C(500)-1.)44,8,44 

IF(C(500)-2.)7,6,7 

C(404)=0. 

Y{231)=0. 

GO  TO  5 
C(404)=l. 

Y(231)=l. 

GO  TO  5 


C(404)=-.C(I2) 

TEMPa(-C(12))*Y(212) 

C(I2)=(Y(211)*I'EXP(TEMP))+(Y(213)1‘C(I2)) 

C(405)=C(I2) 

IF(C(404))45. 19,45 

TEMP=((Y(207)-Y(222))*{y(207)-Y{222)))+4.*Y(225) 
Y(12)=(((Y{207)-Y(222))*(Y(207)-Y(222)))+2.*Y(225))/SQR(TEMP) 
TEMP=(Y(222)*  Y{222))+4.  ♦Y(225) 

TEMPa((Y(i22]*Y(222))+2.»«Y(  225)) /SQR(  TEMP) 
Y{I2)=(Y(207)+(2.*Y(208))-TEMP-Y(I2))/(2.*Y{208)) 
TEMP=C(404)/(((Y(222)-Y{210))*(Y(222)-Y(210)))+C(402)) 
C(12)=(2.  ♦y(208)/(l.  -Y(I2)*C(I2)))*SQR(TEMP) 

IF(SENSE  SWITCH  1)46,  19 
PUNCH,  Y(I2),C(I2),I2,C(405) 

IF(C(404)-1.)47,  11,47 

Y(I2)=Y(225)/(((Y(207)-Y(2Z2))*(Y(207).Y{222)))+4.*Y(225)) 

V(I2)=l.+Y(I2)+(2.25+Y(I2)+Y{I2)) 

TEMP={({Y{207)-Y{222))*(Y{207)-Y{222)))+4.*Y(225)) 

Y(I2)=Y(I2)*(Y(207)-Y(222))/SQR(TEMP) 

C(406)=Y(I2) 

Y(I2)«y(225)/((Y{222)*Y(222))+4.*Y(225)) 

TEMP=<  Y(222)*Y(222))44.  *Y(225) 


Y(I2)=(  1.  +Y{I2)+(  2.  25*Y(I2)*Y(I2)))*(  Y(222)  /SQR(TEMP)) 

IF(SENSE  SWITCH  1)48,49 
PUNCH,  Y(I2),C(406) 

Y(I2)=i(Y(I2)+C(406))*0.  5 

IF(C(404))26,26,27 

C(I2)=0. 

C(I2)=C(I2)+(1.  -SQR(C<404)))/(1.  -Y(i2)*C(40S)) 

IF(SENSE  SWITCH*! )50,  11 
PUNCH,  Y{I2),C(I2) 

TEMP=(((Y(222)-Y(210))*(Y(222)-Y(210)))+C(402)) 
C(I2)=(C(I2)>*'Y(206)*C(405)*C(402))/((C(403)*C(403)*C(403))*TEMP) 
IF(SENSE  SWITCH  1)51,52 
PUNCH,  Y{I2),C(I2),C(404),C{402) 

12=12+1 


Y{218)=Y{218)+Y(217) 

CONTINUE 

Y(219)=0. 

M6-11  - 1 

DO  4  12=3,  M6,  2 
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J»I2+1 

Y(219)=Y(219)+(2.*C(I2))+4.*C(J) 

4  CONTINUE 

J 1=15+500 
J2=I1+1 

Y(Jl)=(Y(2J7)/3.HC(l)+(4.*C(2))+y(219)+C(j2)) 

Y(221)v(y(205)-Y(204))/D4 

C(502)=Y(502) 

C{J1)=Y(J1) 

IF{SENSE  SWITCH  1)53,54 

53  PUNCH,  C{J1),I5 

54  Y(222)nY(222)+Y(221) 

23  CONTINUE 
Y(223)n0. 

DO  24  15=4,14,2 

J4rsI5+500 

J5=I5+501 

24  Y{223)=Y<223)+(2.+C(J4})+4.*C(J5) 

IF(SENSE  SWITCH  1)55,56 

55  PUNCH,  Y(223) 

56  j6«14+2 

Y(224U(  i.(221)/3.  )*(C{502)+(4.  ♦C(503))+Y(223)+C( J6)) 
Y(233)=Y/224)+Y(232) 

PUNCH.  Y(233),Y(204),Y(205),Y(224) 

PUNCH,  y(207),Y(208},Y(209),Y(210) 

IF(SENSE  SWITCH  1)57,58 

57  PUNCH,  Y(229),Y(230),Y(231),C|409) 

PUNCHi  Y{217),  y(221),  y(502),  C(502) 

58  PUNCH,  Y(211),Y(212),Y(213),Y{206) 

J7=0 

PUNCH,  J7 
GO  TO  1 
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